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Abstract 

In this paper, we explore the deep inelastic structure functions of hadrons 
nonperturbatively in an inverse power expansion of the light-front energy of 
the probe in the framework of light-front QCD. We arrive at the general ex- 
pressions for various structure functions as the Fourier transform of matrix 
elements of different components of bilocal vector and axial vector currents 
on the light-front in a straightforward manner. The complexities of the struc- 
ture functions are mainly carried by the multi-parton wave functions of the 
hadrons, while, the bilocal currents have a dynamically dependent yet simple 
structure on the light-front in this description. We also present a novel analy- 
sis of the power corrections based on light-front power counting which resolves 
some ambiguities of the conventional twist analysis in deep inelastic processes. 
Further, the factorization theorem and the scale evolution of the structure 
functions are presented in this formalism by using old-fashioned light-front 
time-ordered perturbation theory with multi-parton wave functions. Nonper- 
turbative QCD dynamics underlying the structure functions can be explored 
in the same framework. Once the nonperturbative multi-parton wave func- 
tions are known from low-energy light-front QCD, a complete description of 
deep inelastic structure functions can be realized. 



PACS numbers: 12.38.-t, 13.85.Hd, 13.88.-Fe, ll.30.Rd 
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I. INTRODUCTION 



With the recent and planned experiments on polarized and unpolarized structure func- 
tions, the field of deep inelastic scattering has entered a new era. New experiments are 
beginning to provide invaluable information on the so-called "higher twist" (power sup- 
pressed) contributions to deep inelastic cross sections, a theoretical understanding of which 
requires nonpcrturbativc information on the structure of hadrons. To unravel this structure, 
there is an urgent need to develop nonperturbative theoretical tools which are preferably 
based on physical intuition and which at the same time employs well-defined field theoreti- 
cal calculational procedures. Towards this goal, in this work, we propose a new method of 
calculation of structure functions combining the coordinate space approach based on light- 
front current algebra techniques and the momentum space approach based on Fock space 
expansion methods in light-front theory in a Hamiltonian QCD framework. 

To get an intuitive picture of deep inelastic scattering in field theory, it is extremely help- 
ful to keep close contact with parton ideas. However, partons were originally introduced as 
coUinear, massless, on-mass shell objects. In reality, the QCD governed interacting partons 
should not be collinear and massless. The question, then, is can one generalize this concept 
and introduce field theoretic partons as non-collinear and massive (in the case of quarks) but 
still on-mass shell objects in interacting field theory? Nonperturbative hght-front Hamil- 
tonian description of composite systems which utilizes many body wave functions for the 
constituents allows us to precisely achieve this goal. 

Now, the problem is how to introduce many body (or multi-parton) wave functions in the 
description of deep inelastic structure functions. An attractive possible avenue is provided 
by the BJL (high energy) expansion of scattering amplitudes together with the use of light- 
front current algebra. This is essentially a nonperturbative approach where the expansion 
parameter is the inverse of the light-front energy of the probe (in the present case, the virtual 
photon). In this approach one can arrive at expressions for various structure functions as the 
Fourier Transform (in the light-front longitudinal direction) of matrix elements of different 
components of bifocal vector and axial vector currents in light-front field theory. 

In the standard approach to deep inelastic scattering based on Wilson's Operator Prod- 
uct Expansion (OPE) method which is more mathematical, one considers the problem of 
renormalizing composite operators. In contrast, in the current algebra approach, products 
of field operators are evaluated at equal light-front time and have the same form as their 
free field theory counterparts. Hence in this approach, most of the complexities appear to 
be carried by the hadronic states. We demonstrate that an expansion of the state allows 
us to exhibit this complexity manifestly in terms of the multi-parton wave functions, where 
the constituents are on- mass shell objects with non- vanishing transverse momenta. The 
nonperturbative nature of the matrix elements relevant for various structure functions are 
thus directly translated into the language of multi-parton wave functions. The structure 
functions, then, can be conveniently evaluated once the nonperturbative wave functions 
renormalized at scale Q are known. At present major efforts are underway to evaluate these 
wave functions. 

With the advent of QCD, the current algebra approach which originally lead to the pre- 
diction of scaling without introducing the concept of partons was altogether abandoned in 
favor of the OPE method, primarily because it was known that the canonical manipulations 
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that lead to many of the current algebra predictions are invalidated in perturbation theory 
which gives rise to logarithmic corrections. The only exceptions were sum rules protected 
by conservation laws. Obviously, what was missing in the current algebra approach was the 
realization that the matrix elements relevant for deep inelastic scattering are those renormal- 
ized at the physical scale Q. Thus one of our major problem is whether we can consistently 
carry out renormalization procedure. In a subsequent paper we shall demonstrate that 
in the case of leading order structure functions the scaling violations of perturbation theory 
can be successfully addressed by replacing the hadron target by a dressed parton target in 
the matrix element and carrying out a well-defined perturbative expansion which closely 
follows the techniques of (light-front) time-ordered perturbation theory [0-^. Here, we will 
show that we can address the issues of factorization and scale evolution in the case of a 
hadron target by separating the soft and hard parts of the multi-parton wave functions. 
Furthermore, the nonperturbative contributions to deep inelastic structure functions can 
be addressed within the same framework by incorporating the newly developed light-front 
renormalization group approach to nonperturbative QCD [^-0]. Therefore, a unified light- 
front description of the perturbative and nonperturbative QCD underlying the deep inelastic 
structure functions can be realized. 

The paper is organized as follows. In Sec. II, a brief overview of deep inelastic structure 
functions is presented. In Sec. Ill, we discuss within light-front QCD how the deep inelastic 
structure functions are nonperturbatively related to the bilocal operators in light-front cur- 
rent algebra through an inverse power expansion of light-front energy of the probe. In Sec. 
IV, we determine the structure functions in terms of the matrix elements of the light-front 
bilocal vector and axial vector currents separated only in the longitudinal direction. The 
operator structures involved in the structure functions are explored and the complexities of 
structure functions are analyzed in Sec. V, where a scheme to evaluate the hard and soft 
contributions to deep inelastic structure functions in the light-front time-ordering Hamilto- 
nian formalism is proposed. Meanwhile, the concept of twist is examined and redefined on 
the light-front, which removes some ambiguities in previous works. In Sec. VI, we re-derive 
sum rules the structure functions obey and discuss their physical implications. Finally a 
summary is given in Sec. VII. 

II. A BRIEF OVERVIEW ON DEEP INELASTIC STRUCTURE FUNCTIONS 

We begin with a brief review of the basic ingredients of lepton-nucleon deep inelastic 
scattering (DIS): 

e{k) + h{P) —^e{k-q)+X{P + q). (2.1) 
The cross section for the above scattering process is given by 

da 1 a^E' 



dVtdE' 2M g4 E 



L,.W^\ (2.2) 



where E [E') is the energy of the incoming (outgoing) lepton, L^j^^ is the leptonic matrix 
element. 
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= ^ ^hM^'^ s')u{k', s')^^u{k, s)] 

s' 

— 2{k kn + kj^kfj) — 2gi^i,k ■ k — 2ief^i,p^q'^s'^ , 



(2.3) 



and W^'^ the hadronic tensor which contains all the hadronic dynamics involved in DIS 
process, 



An 



d'^ e'^<{ps\[j^{0,rm\PS), 



(2.4) 



where P and S are the target four-momentum and polarization vector respectively (P^ = 
M'^,S'^ — —M'^,S • P — 0), q is the virtual-photon four momentum, and J'^{x) — 
Ha^a'^ai.^)l^''Pa{.x) the elcctromagnetic current with quark field ipai^) carrying the flavor 
index a and the charge e^. 

The above hadronic tensor can be decomposed into independent Lorentz invariant func- 
tions: 



-ie^^^^qx 



S^Ws{x,Q^)+P,S-qW^{x,Q^) 
I (g'" - ^) F,(x, Q') + \P^P^ - 4 (P^q^ + P^q^) + g^^"^] ^^^^ 



V 



SaLgi{x, Q^) + SaTgrix, Q^) 



The dimensionless functions 



Fl{x,Q^)^2 -Wi + [M^- 



{P.qf 



(2.5) 



(2.6) 



and 

F2{x,Q'')^vW2{x,Q'') (2.7) 
are the so-called unpolarized structure functions measured from the unpolarized target and 



and 



yi(x, g") = V w^{x, g") + vWi{x, g") 



grix, g^) = g^{x, Q^) + g^ix, Q^) = vW^{x, g^) , 



(2.8) 
(2.9) 

are the longitudinal and transverse polarized structure functions. While, a; = ^ is the 
famous Bjorken scaling variable, — —q^ the momentum transfer carried by the virtual 
photon and v — P ■ q the energy transfer. The longitudinal and transverse polarization 
vector component are given by 



Sjx s. 



s. 



HT — Sfj, — P/j, 



S-q 



V 



(2.10) 



These structure functions provide a probe to explore various aspects of the intrinsic 
structure of hadrons. It may be worth noting that in the hterature, gi and g2 are usually 
used to characterize the longitudinal and transverse polarized structure functions. However, 
g2 is not really a transverse polarized structure function. It also has no clear physical 
interpretation. Only gx which can be directly measured when the target is polarized along 
the transverse direction characterizes the full information on the transverse polarization 
structure. 

The early SLAG experiment discovered that the structure function F2 depends only 
on the Bjorken scaling variable x, and is independent of momentum transfer Q^. This 
discovery lead to the parton picture proposed first by Feynman ^ in which the constituents 
in hadrons can be treated as point-like free particles, i.e., the partons. Of course, the picture 
of point-like and non-interacting partons is certainly over-simplified. As it is well known now, 
the fundamental constituents (i.e., quarks and gluons) inside hadrons are indeed strongly 
interacting with each other and are governed by the fundamental QCD theory. Only in the 
limit of very large Q^, the asymptotic freedom feature of QCD makes the quarks and gluons 
behave as point-like and weakly interacting partons. The scale invariance of F2 is violated 
for finite values of Q^. For the past twenty- years, QCD investigations on deep inelastic 
structure functions are mainly concentrated on the behavior of the scaling violation of these 
structure functions, which is believed to be the best way to test QCD as the fundamental 
theory of the strong interaction. 

Among these investigations, two methods have dominated the whole research topic: the 
OPE method and the QCD improved (or field theoretic) parton model based on the 
factorization scheme [|T0|. The structure functions have been extensively explored either 



in terms of its moments in the language of OPE or in terms of diagrammatic calculation 
in the language of QCD improved parton model. Both the methods have made great suc- 
cess in understanding the scale evolution of the structure functions within the perturbative 
QCD domain. Although the deep inelastic scatterings provide a novel way to separate the 
physics of perturbative and nonperturbative QCD dynamics and allow the exploration of 
high energy behavior of the constituents inside the hadrons, the structure functions them- 
selves, however, are still dominated by the nonperturbative strong interaction dynamics. A 
complete understanding of hadrons crucially depends on our understanding of nonperturba- 
tive QCD. The OPE method addresses the structure functions in terms of their moments, 
which naturally separates the short- distance and long-distance dynamics but it also makes 
the nonperturbative dynamics of long-distance physics more complicated in terms of the 
moments. On the other hand, the QCD improved parton model was built in the framework 
of perturbation theory with the assumption of coUinearity, which simplifies the perturbative 
QCD treatment but it is also unclear how to explore nonperturbative QCD dynamics be- 
cause the nonperturbative QCD dynamics is mainly determined by the non-coUinear motion 
of the low energy quarks and gluons. 

In the next section we shall use an inverse power expansion of the light-front energy 
q~ of the virtual photon (the extended BJL theorem on the light-front) to extract these 
deep inelastic structure functions. This approach was originally proposed to study DIS 
sum rules protected by conservation laws in the pre-QCD era [|11[]. Here, we shall extend 
this approach to QCD without recourse to perturbation theory. Therefore it is essentially 
a direct nonperturbative QCD description. The resulting structure functions are directly 



5 



expressed in terms of the hadronic matrix elements of light-front bilocal vector and axial 
vector currents, where the currents have a relatively simple structure although they are 
dynamically dependent. All the hadron dynamics (including both the perturbative and 
nonperturbative dynamics) reside in the multi-parton hadronic wave functions. This is a 
description very different from that of the OPE method and the QCD improved parton 
model. 



III. AN EXPANSION IN INVERSE POWER OF LIGHT-FRONT ENERGY OF 

THE VIRTUAL PHOTON 

The inverse power expansion of the virtual photon light-front energy q~ is applied to 
forward scattering amplitudes. Explicitly, as it is well known, the hadronic tensor is related 
to the forward virtual-photon hadron Compton scattering amplitude: 

W^"" = —ImT^" (3.1) 
27r 



with 



S^T,{x,Q') + P^S-qT,{x,Q') . (3.2) 



Using the optical theorem, we have 



Ux, Q') = 2 r dq'^ ^f/'^J , ^ = 1, 2, 3, 4 . (3.3) 

J— oo q — q 

Then, as we shall see, the structure functions are connected with the light-front bilocal 
currents through the l/q~ expansion of T'^'^. 



An expansion oiT^^ in terms of 1/g was originally proposed by Jackiw et al. ||TT| based 
on BJL theorem. The general expansion in 1/q^ is given by 



oo -1 I 1 /• 

^'^'^ = - E (— )" / dCd''Ue''^<{PS\[[id^T.P{i), rmi^=o\PS) , (3.4) 



1 \n+l 

n=0 

where q~ = q^—q^, the light-front energy of the virtual photon, and d~ = 2^^ is a light- front 
time derivative. The light-front coordinates of the space-time are defined by 



e = e±e , ei = e(^ = i,2). (3.5) 

The above expansion shows that the time-ordered matrix element can be expanded in terms 
of an infinite series of equal light-front time commutators. 

For large and large v limits in DIS, theoretically without loss of generality we can 
always select a Lorentz frame such that the light-front energy q~ of the virtual photon 
becomes very large. Then, only the leading term in the above expansion is dominant, i.e.. 
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J^flU 



largo q 



— I drd'Ue''-^{PS\[J^{0,J''mi^=o\PS). 



(3.6) 



As a result, the leading contribution to the deep inelastic structure functions is determined by 
the light-front current algebra. The light-front current commutator can be computed directly 
and exactly from QCD (where QCD should be quantized on the light-front time surface = 
^0 + ^3 _ g ^j^]^ light-front gauge ^4+ = [0,0). Hence, all the subsequent derivations 
are exact within the light-front QCD and without further assumptions or approximations of 
the collinear and massless partons that were used in the previous derivations [^,10 . 
Explicitly, the basic commutation relation on light-front is 



{^+(x) , i>i{y)}a,+=y+ = A+(5(x - y )(5^(xj_ - y±) 



(3.7) 



which is exact in the full QCD theory [Q, where ip+{x) is called the dynamical component 
of fermion field on the light-front: 



1 



^(a;)=V^+(x)+^_(x) , ^±{x)=A^^{x) , A± = ^7 V • (3-8) 
The minus component ip^lx) is determined from ?/'+( result of Dirac equation: 

ip-{x) = -^{ia± ■ + I3m^ilj+{x) . (3.9) 

d± — igA± is the 



Here we have already used the light-front gauge = 0, and D± 
transverse component of the covariant derivative, = 7^7*, (3 = 7*^. 

Because of the above special property of quark (or more generally fermion) field on the 
light-front, the light-front current explicitly depends on interaction of the theory, which 
is very different from the usual equal-time formulation. In other words, the fundamental 
interaction is manifested explicitly in the light-front current commutators. 

From Eqs.( p.7|) and ( p.9|) , we have 



{^+(x) , r-{y)}. 



A+ 

17 



■e{x - y 



ia±_ ■ D]_ — I3m 5'^{x<^ — y^) 



(3.10) 



Thus, after a tedious but straightforward calculation, one can find that 



J+(x) , J-{y) 



1 



{x -y )S^{x± - y±)V^ {x\ 



+di [-e{x- - y-)6'{x^ - y^) [vi{x\y) 



-^le'^A^^y) 



— h.c. 



(3.11) 



where and are defined as the bilocal vector and axial vector currents: 

v^ix\y)=^^ix)rMy), 



(3.12) 
(3.13) 
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As we can see, the light-front current commutators are very different from the equal- 
time current commutators. Here the commutator is indeed given by terms involving spatial 
derivatives. These space-derivatives come from the non-locality of ip-{x) on the light-front. 
In the equal-time formulation, there is no such nonlocality to the fermion field. Therefore 
one cannot derive such a commutator from the naive canonical equal-time commutations. 
As we will soon see in the next section it is these spatial derivatives that lead to the simple 
expressions of the structure functions in terms of bilocal current matrix elements. This is an 
essential feature in the present approach that make the light-front current algebra specially 
useful in the exploration of the deep inelastic structure functions. 

The commutators for other current components can also be found straightforwardly. For 
example. 



J+(x) , J'{y) 



-e(x -y )6\x^-y^)V^{x\ 



a { 

Ux- - y~)5\x^ - y^) {g'^V^{x\y) + ^e^M+(x|y)) 



h.c. . 



(3.14) 



Thus, one can use Eq.( |3.11| ) to extract the structure functions and then use Eq.( p.l4D to 
make a consistency check. 

Now, the Compton scattering amplitude in the large q~ limit can be immediately ex- 
pressed in terms of the hadronic matrix elements of the bilocal vector and axial vector 
currents. For example, the (H — ) component T+~ is given by 



large q 



mps\Y.^i{lq''v-{i-\Q) 



- \qimc\^)+ie''K{rm] - h.c. \ps) . 



(3.15) 



The above result shows that the bilocal vector and axial vector currents entering in T'^'^ are 
separated only in the longitudinal direction ^~ = This property naturally leads to the 

well-known scaling behavior of the structure functions when we ignore the QCD dynamics 
at very high Q^. In the next section, we will extract explicitly the structure functions in 
terms of these bilocal current matrix elements. 



IV. THE GENERALIZED EXPRESSIONS FOR DEEP INELASTIC STRUCTURE 

FUNCTIONS 



Now, let us pick up the same (H — ) component of the hadronic tensor Eq. ( p.5| ) to find 
the deep inelastic structure function by comparing with Eq. ( p.l5| ) through Eq. (|3.3|) . In the 
large q~ limit. 



1 



2 V q'^ 



2le''q^ 



c 9i , o 9t 

DjL h OjT 



(4.1) 



where SjT = Sj — S^-^ and SjL = Sj — Sjt = S^^, and u = \P^q in the large q limit. 
We introduce the form factors for the bilocal current matrix elements. 
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{psiv^m - v^m\ps) = p^v,^{p\^- p) + el"2.(p^e ■ p) , (4.2) 
{ps\A>i^{m + A^Mi)\ps) = s'^A^iP', ^.p)+p^^. sAaiP', e ■ P) 

+es-a3a{p',^-p). (4.3) 

Since, ^^'-^ = 0, it follows that the matrix elements of the plus and transverse components 
of the bilocal current yield the same form factor Via- Using the definition 



e(r) 



TC J~oo g"" 



we find that 



nq 



oo q'~^ — g+ 



e 



-XP-q^ -P±-q±)V 



la 



+ ^q^^V2a + -e''q^[S,A,^ + PJ ^ 



■A 



2a 



The bilocal current form factors are determined from Eg. ([4 .21): 



Vic 



V 



1 

pi 
1 



2a 



Ale 



-(P5iv^„(r)7Va(o) - i^Mfi^a{n\ps) , 

p- 
p+ 
p- 



-{ps\^^{r)[i 



-l^)MO)-h.c.\PS) 
-^')MO)-h.c.\PS), 



S: 



A 



2a 



P+^-S't 



Comparing with Eqs.( |4.1D and (|4.5| ) through Eq. (|3.3|) , we obtain, with t] = |P^^ 

' dve-''"" E e« (^^1 iCh^i^a (0) - (o)7+V^. (r ) I 

dr]e^''^''j:el{PS\^^{nilMO)-^a{OhlMn\PS) , 



47^P^ 
1 



47rP! 



(4.4) 



(4.5) 

(4.6) 
(4.7) 

(4.8) 
(4.9) 
(4.10) 

(4.11) 
(4.12) 
(4.13) 



where the last equality is found for the first time here |]T2[. Its physical interpretation will 
be given later. 

/ dve-^^^ Y: el [(P- - ^)Fi„ + rV^a 



nP+q- 



X 7 



_P1 

:p+ 



7i^a(0)-/l.C.|P5) 



(4.14) 
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where the first equahty may be reduced to the same expression obtained by the collinear 
expansion in the Feynman diagrammatic method up to the order twist-four [|l^ . But here it 
is directly obtained in the leading order in the l/q~ expansion without involving the concept 
of twist expansion. Moreover, the polarized structure functions can also be found directly 

as 

9i{x,Q') = ^J dve-'^^Y.^l{^i» + lp^rA2a) (4.15) 
= ^ / t/r^e-''^Ee'(^^lV^.(r)7+75V^a(0) + V^(0)7+75V^(r)|i^5), (4.16) 
9t{x, Q') = ^ [ rf^e-^'^^ E elA,^ (4.17) 

= 8^/ dve-''"'T.el{PS\i^^ir){Y - ^l^)lnMO) + h.c.\PS) . (4.18) 

The above results are derived without recourse to perturbation theory, and also without 
the use of concept of collinear and massless partons. They are also the most general expres- 
sions for the leading contribution (in the l/g~ expansion, not the leading contribution in 
terms of twists) to the deep inelastic structure functions in which the target is in an arbitrary 
Lorentz frame. Some of these expressions have not ever been obtained in previous works. 
The x-dependence of these structure functions is obvious in the above expressions. The 
scale ((5^)-dependence is hidden in the hadronic bound states \PS) which can be described 
by multi-parton wave functions. In the next section, we shall analyze the complexities of 
structure functions in terms of bound states in our description. 



V. COMPLEXITIES OF STRUCTURE FUNCTIONS 
A. Multi-parton wave functions 

As we have seen, the derivation of structure functions in the previous section is apparently 
quite different from the QCD improved (or field theoretical) parton model. The latter which 
is based on the colhnear concept is purely an application of perturbation theory where 
further exploration of nonperturbative dynamics is lacking. Our description is also obviously 
very different from the OPE method. On the other hand, structure functions themselves 
are dominated by the nonperturbative quark-gluon dynamics. When we formulate them 
on the light-front, the structure functions are proportional to the simple hadronic matrix 
elements of the bifocal currents that are separated only in the longitudinal direction. In this 
formulation, no time evolution or propagation is explicitly involved in the matrix elements. 
Hence, unlike the OPE or the perturbative field theory descriptions of parton model, all the 
perturbative and nonperturbative dynamics here are completely carried by the structure of 
target's bound state. This is closer to the real physical picture probed in experiments. 

The bound state of a hadron on light-front can be simply expanded in terms of the Fock 
states, 

\PS) = J2 J dxi(fK±i\n,XiP^,XiP± + K±i,Xi)^^{xi,K±i,Xi), (5.1) 

n,\i 
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where n represents n constituents contained in the Fock state |n, XiP~^, XiP± + Aj), Aj is 
the hehcity of the i-th constituent, /' denotes the integral over the space: 



^ = 1, and 'YK±i = 



(5.2) 



while Xi is the fraction of the total longitudinal momentum carried by the i-th constituent, 
and K±i is its relative transverse momentum with respect to the center mass frame: 



Xi 



P+ 



Ki± = Pi± - XiP± 



(5.3) 



with pf, pi± the longitudinal and transverse momenta of the i-th constituent. $f (xj, hi±i, Xi) 
is the amplitude of the Fock state \n,XiP^ ,XiP± + K_Li,Ai), i.e., the multi-parton wave 
Junction, which is boost invariant and satisfy the normalization condition: 



dXid^K^i\^l{Xi,K^i,\i)\^ = 1, 

n,Xi 



(5.4) 



and is, in principle, determined from the light-front bound state equation. 



i=l 



Q1Q9 



{qqq\Hint\qqq) {qqq\Hint\qqqg) 
{qqqg\Hint\qqq) 



$5 
qqq 



(5.5) 



Here Hint is the interaction part of the light-front QCD Hamiltonian 0]. Thus, the complex- 
ities of the structure functions carried by hadronic bound states are now translated into the 
language of multi-parton wave functions on the light-front, rather than composite operators 
in OPE. 

Explicitly, let us look at the structure function F2{x, Q^). It is found for the first time by 
us 1 12 1 that F2 can be expressed in terms of a matrix element of either the plus component 



or the transverse component of the bifocal vector current. On the light-front, these two 
components have totally different operator structures but amazingly their matrix elements 
determine the same structure function. 

The plus component (usually called the "good" component), 



(5.6) 



has no explicit dynamical dependence, and has the lowest mass dimension (a twist-two 
operator in OPE language). The corresponding matrix element has straightforward parton 
interpretation. It is clear from the above operator that on the light-front it is just a quark 
(parton) number operator which immediately leads to the fact that F2 is proportional to 
parton density distributions qa{x,Q^). 



^^^^^ = Y.^lq.M'). 

X Q 

g«(x,Q') = / d^k^{PS\Y,hi{k,x)h^{k,\)\PS) 

= j d^K^Yl j dXid^Ku\^n^^{x,Xi,KuAi)\' 
n.Xi 



(5.7) 



(5.8) 
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where the Q^-dependence is carried by the multi-parton wave functions with the active 
parton renormahzed at the scale Q^, /" denotes the integral in the right-hand-side over the 
space of Eq. ( |5^ ) except for the active parton {x,k±) = {k~^ /P~^ ,p± — xP±). With this 
consideration it is straightforward to derive the logarithmic corrections that is the same as 
that obtained in the QCD improved parton model or in the OPE, as will be given in . In 
this case, all the three descriptions are almost the same. The only difference here is that in 
our framework, the perturbative QCD dynamics is transferred from the composite operator 
into the scale-dependent multi-parton wave functions on the light-front, which enables us to 
describe the nonperturbative dynamics in the same framework. 

But the transverse component (sometimes called the "bad" component), 

= V^-7iV^+ + ^+71 V'- , (5.9) 

depends explicitly on the fundamental quark-gluon interaction in QCD. According to the 
twist analysis [|13|, the transverse component of the bilocal current is a twist-three operator 



which has no simple parton interpretation. However, we have explicitly shown that 
the corresponding matrix element of the above transverse component must have the same 
parton interpretation as that from the plus component. This is indeed obvious because 
they represent the same form factor of the bilocal current [see Eqs.( |4.6| ) and ( [4.7| )] and they 
describe the same structure function F2 [see Eqs. (|4.12| ) and (|4.13|) ]. This conclusion is very 



different from the current understanding of the twist of composite operators in the OPE or 
in the QCD improved parton model. 



The explicit calculations in Ref. |]12| further demonstrate that the real dynamics con- 
tained in the structure functions is determined by the matrix element with the rich in- 
formation carried by the multi-parton wave functions. It is the complicated multi-parton 
wave functions that causes the same behavior for the matrix elements of the two apparently 
different operators (the plus and transverse components of bilocal current). However, this 
property was overlooked in previous works. Partly because one usually worked in a spe- 
cific Lorentz frame, such as the rest frame or the infinite momentum frame in which the 
target's transverse momentum P± = 0, and also because enough attention was not paid on 
the hadronic bound state structure. With P± = 0, F2 can only be expressed by Eq. ( [4.121) . 
The expression in terms of the matrix element of the transverse component of the bilocal 
current, i.e. Eq. ( [4.13| ), does not exist if one lets Pi = 0. In other words, picking up a 



specific Lorentz frame on light-front may lead to some ambiguities. In the next subsection, 
we shall present some discussion on such ambiguity in the concept of twist which is currently 
one of the most interesting topics in the study of deep inelastic processes. 



B. A Lorentz invariant definition of the twist 



As we have seen the lack of the additional expression ( [4. 13] ) for F2 in other approaches 
may be due to the specific choice of the Lorentz frame one often used. Indeed, since the deep 
inelastic processes are dominated by the physics close to the light-front, one often choose 



the following parameterization to analyze the structure functions [10,13 



(5.10) 
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where p'^ = (p, 0,0,p) and = (2,0,0, —2) are the hght-hke vectors; p = M/2 yields the 
target rest frame and p ^ oo leads to the infinite momentum frame. This choice of Lorentz 
frame for the deep inelastic processes may simplify the analysis, but it loses the generality 
when one tries to extract some general conclusions. A typical example is the twist analysis 
in the determination of the dominant contributions to deep inelastic structure functions. 

The dominant contribution to the structure functions in terms of twist expansion is 
the basis of the OPE and QCD improved parton model analysis. Measuring the higher- 
twist contributions has become a main topic in the current and future experiments since it 
may provide some nontrivial test of perturbative QCD beyond the logarithmic corrections. 
However, the definition of twist in literature is not unique and its physical interpretation 
is also not very transparent. According to standard definition, twist of a given composite 
operator is given by mass dimension of the operator minus its spin. In reality, this definition 
has no direct connection with experimental parameters. 

Practically, the twist can be given in a less formal way: the twist of an invariant matrix 
element of a light-front bilocal operator, which controls the scaling behavior of the matrix 
element is determined by the power of l/Q that it contributes to deep inelastic processes 
||T3| . This definition is correct if the invariant matrix elements here are specified by those 
Lorentz invariant form factors of the bilocal currents given in Eqs.( |4.6[ - [4.m| ). However, in 
previous applications, one used the specific coordinates ( p. 10 ) to define the invariant matrix 
elements of the light-front bilocal currents and further used the usual dimensional analysis to 
determine the corresponding twists. Unfortunately, the twist of an matrix element defined 
in such specific coordinates with the usual dimensional analysis is not Lorentz invariant and 
contains ambiguous. 

For example, according to the above definitions, Jaffe and Ji concluded that plus com- 
ponent of bilocal current is a twist-two operator, and transverse and minus components are 
twist three and four operators [|l^. As a result, matrix elements of the transverse component 
of bilocal currents should be suppressed by a power of 1/Q in comparison to the plus com- 
ponent matrix element. However, we find that the transverse component matrix element is 
actually suppressed by P± rather than Q^: Smaller the value P± takes in a specific choice of 
a Lorentz frame for the target, smaller will be the contribution arising from the transverse 
component matrix element. In the infinite momentum frame or in the target rest frame 
where the target has zero P±, the transverse component matrix element must vanish. On 
the other hand, the plus component matrix element retains its value, no matter what the 

value is (modulo logarithms). Thus, one cannot obtain the conclusion that transverse 
component of bilocal current is a twist-three operator that is suppressed by the factor 1/Q 
in hard processes. The fact that the transverse component matrix element vanishes is a 
consequence of using the specific Lorentz frames with P± = and is independent of the 
value of Q^. 

Indeed, we also find that twist of Lorentz invariant matrix elements of light-front bilocal 
currents can be properly defined in terms of light-front power counting analysis for light-cone 
dominated deep inelastic processes. Here one needs to consider the target in an arbitrary 
Lorentz frame. Then Eqs.( [4.2|) -( [4.3|) give the general definition of the Lorentz invariant 
matrix elements, i.e., the form factors of the bilocal vector and axial vector currents. The 
twist of a Lorentz invariant matrix element of bilocal currents is thus directly given by 
the mass dimension of the corresponding bilocal current form factors in light-front power 
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counting. One must notice that light-front power counting is very different from the usual 
one 1^. The hght-front power counting separately introduces different scales for transverse 
and longitudinal dimensions. Only the transverse dimension has mass scale, 

m : — . (5.11) 

x± 

In light-front power counting both the light-front quark field 4'+ and the gluon field A± 
have mass dimension one because their power assignments are given by 

: -L— , : — , (5.12) 

although their longitudinal dimensions are different. Besides, the transverse derivative has 
mass dimension one and the longitudinal derivative has no mass dimension. Thus, from Eq. 
( p.9|) , we can see that ip^ : has mass dimension 2. 

In deep inelastic scattering processes, we have shown in Eqs. (|4.2|) and (|4.3|) that there 
are five Lorentz invariant form factors associated with the matrix elements of the bifocal 
vector and axial vector currents: Via, V2a, ^la, and A^a- With simple power counting, 
one can easily find that three of them. Via, and A2a, are twist-two, while the remaining 
two, A2a and A^a, are twist-four. Based on the normalization of bound state 

{P'S'\PS) = 2(27r)3p+(5'^(P - P')6ss' , (5.13) 

the light-front mass dimension of the state \PS) is —1. Then we can easily count the mass 
dimension for each bifocal current form factor in light-front power counting as follows: 



Vi_ 


~ {PS\ 




^+ 


\PS) 




m-dim(\/i) 


= -1 


+1 


+1 


-1 
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m-dim(A2) = 
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{PS\ ^1 tp. \PS) 


m-dim(\/2) 
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A,_ 




{PS\ \PS) 


m-dim(y43) 




-1 +2+2-1=2. 



The twist of the Lorentz invariant matrix elements (i.e. the bilocal current form factors) is 
defined as its mass dimension plus 2. Therefore, Via, Aia and A2a have twist-two, and V2a 
and A^a are twist-four Lorentz invariant matrix elements. 
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As a matter of fact, there is indeed no twist-three Lorentz invariant matrix element in the 
leading contributions to deep-inelastic lepton-hadron processes. This conclusion should not 
be very surprising since experimentally no twist-three matrix element with a 1/Q coefficient 
has been extracted in deep inelastic processes. Theoretically, it is also hard to find some 
result with power correction proportional tol/Q rather than l/Q^ in inclusive lepton-hadron 
deep inelastic processes. Our analysis shows that the next higher twist corrections to the 
leading twist-two contributions in deep inelastic processes are twist-four contributions which 
should be suppressed by the factor l/Q^. Indeed, from light-front power counting analysis, 
the mass dimension of Lorentz invariant bilocal current form factors is much more meaningful 
than the conventional definition of twist. If the light-front mass dimension of a form factor 
is zero, it means that there is no l/Q power suppression, see Eqs.( [4.liD , ( [4.15| ), and |4.17|) for 
F2, gi and qt- If the mass dimension of a form factor is not zero, then there might be some 
0{1/Q) power suppression to ensure that the structure function has the same dimensions 
[cf. Eq.(|mi) for Fl]. _ 

Now we can see that because the Lorentz invariant form factor Via is a twist-two matrix 
element, no matter if it is expressed in terms of the matrix element of plus or transverse 
component of the bilocal vector current, the leading contribution of F2 is always a twist-two 
contribution, as is already known. The contribution of the transverse component matrix 
element goes to zero when we analyze it in the specific light-cone coordinates ( |5.10|) , but 
the ratio of the transverse component of the matrix element and P±, which is proportional 
to F2, remains invariant, as a result of Lorentz invariance. Thus the possible inconsistency 
of eqs. ([4.12|) and ( [4.13|) in terms of the naive definition of twist does not really exist. It 
appears that the conventional definition of twist that is currently used in the literature may 
not be proper in the analysis of deep inelastic processes. 

It is also to be noted that the twist analysis is meaningful only in some restricted portion 
of the phase space. For example, in the region of phase space connected to the limit a; — > 1, 
that is when one of the particle in the Fock state is carrying almost all the hadron momenta, 
the states are far off its energy shell. Apparently, it seems that the contribution from such 
states to the structure functions are power suppressed, for it needs multiparton interaction 
in the first place to dump all the longitudinal momenta to one parton to produce such a 
Fock state. Yet their contributions are comparable to the leading twist ones where center of 
mass energy remains fixed. In fact, these contributions are more important and controls the 
behaviors of structure functions in the x — >■ 1 limit ||T5|Jl6[| . In any case, as shown in these 
references, the twist expansion looses its significance all together in this limit. 

Next, we shall analyze the twist of polarized structure functions, which has remained 
ambiguous for many years in the literature. 



C. Twist analysis of the polarized structure functions 

With the above analysis of the complexities of multi-parton wave functions and the 
new understanding of the concept of twist, we now explore the properties of the polarized 
structure functions gi and gr- 
in the literature, the properties of g2 are usually discussed rather than the real transverse 
polarized structure function g^ for historical reasons. But as we have pointed out, g2 by 
itself has no real physical meaning. However, to clarify the confusion about the contributions 
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of different twists to g2 in the literature, we focus our attention on g2 here. The polarized 
structure functions, especially the transverse polarized structure function, have recently 
received much theoretical and experimental attention. The polarized structure functions 
Qi and g2 are of great interest due to the fact that the former has been considered the 
best quantity to understand the origin of the proton spin, while the later may lead to the 
first observation of higher twist (twist-three) contribution in DIS. In the early naive parton 
model, Feynman claimed that g2, just like gi, has a simple parton interpretation because 
grix) = gi{x) + g2{x) = ^J^q^l^Qrix), where AqT^x) is the distribution of transversely 
polarized quarks in a transversely polarized nucleon Wandzura and Wilczek studied 
the properties of g2 in OPE analysis ||T^ and they claimed that except for a twist-three 



contribution (x) which may be negligible in their model calculation, g2{x) can be related 
to an integral over gi{x): 

92{x) = gf'^{x)+U^), ^(x) = -gr{x) + (5-14) 

Jx y 

The relation between g^'^ and gi{x) is called Wandzura- Wilczek (WW) relation. Later, 
Shuryak and Vainshtein pointed out that the twist-three contribution ^2 is a direct quark- 
gluon interaction effect which is important for g2{x, Q^) so that the measurement of g2 
may also be very sensitive to the interaction dependent higher twist effects in QCD |13| . 

Here we shall apply light-front power counting to the analysis of the twist of polarized 
structure functions. The operator involved in gi is relatively simple, it is the plus component 
of bilocal axial vector current, 

^aiCh^l^i^a = 2^i(r)75^+(0) , (5.15) 

which is a twist-two operator in the conventional definition. The corresponding Lorentz 
invariant matrix element is also twist-two in our definition. However, the operator deter- 
mining the transverse polarized structure function g2 is much more complicated. From Eqs. 
( [4.1(j| ) and ( ^4.181 ), one can obtain: 



92{x,Q')=gT{x,Q')-gi{x,Q') = J c^rye-^^^Ee'^^s. 

^ / rfr/e-''^Eei(^^l[^«(r)(7' - 1^7^)75^^.(0) 



+ h.c. \PS). (5.16) 



The operator in the above matrix element is a mixture of operators of the transverse and 
plus components of bilocal axial vector current and therefore involves twist-two and twist- 
three contributions in the conventional definition. On the light-front, it can be expressed 
as 

?(r)(7± - |t7+)75V^(0) = 1^0^ + 0^ + O'- + (5.17) 

which corresponds to the conventional defining twist-two operator and twist-three op- 
erators of the quark mass part O"^, the quark transverse momentum part O''^ and the 
quark-gluon coupling part O^, 
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O^ = -2^i(r)75^+(0), 

= m4(r)7±(^ - ^)75^+(0) , 

O'- = -ijl(C)h±^f± + ^^7±)75^+(0) , 
a+ (7+ 

o' = 9^i{n(Mn^i±-i±^Mo))i5Mo) (5.18) 

where m and are the quark mass and quark-gluon couphng constant in QCD, and A± = 
A']_Ta is the transverse gauge field. The main concern in the previous works is which operator 
dominates g2 in hard processes in terms of twist analysis. 

However, Eqs. ( |5.16| ) and ( [4.16| ) show that qt is determined by Lorentz invariant form 



factor of the bilocal axial vector current, Aia, while gi is determined by both and 
Ala- According to our analysis on twist, both Aia and 742a are twist-two invariant matrix 
elements because their mass dimensions are zero. Although these operators are apparently 
so complicated, their contributions (in terms of matrix elements) to g2 are indeed twist-two 
regardless if they depend on the quark mass, quark transverse momentum and the direct 
quark-gluon coupling. As a matter of fact, it is not meaningful to separate g2 into twist-two 
and twist-three parts as in Eq. ( p.l4| ) . In our recent work on the validity of Wandzura-Wilczek 
relation, we have shown that such a separation is strongly violated in perturbative QCD. 
The violation originates from the cancellation between the naive twist-two and twist-three 
contributions due to the multi-parton structure of wave functions [|l^. Such a cancellation 



is hardly understood in the conventional twist expansion if one believes that the higher twist 
is suppressed by an order of 1/Q in comparison with the lower twist. Now, we understand 
that there would be no such ambiguity in the twist analysis contributions of g2, if one used 
the definition of twist based on light-front power counting and the Lorentz invariant matrix 
elements of the bilocal currents (rather than the bilocal current operators themselves). 

In conclusion, the complexities of deep inelastic structure functions are mainly carried 
by the multi-parton wave functions of the hadrons which completely determine the Lorentz 
invaraint matrix element (i.e. the form factors) of the bilocal currents. 



D. A scheme for the evaluation of soft and hard contributions to deep inelastic 

structure functions 

Up to this point, all the derivations and discussions of the deep inelastic structure func- 
tions in the l/q~ expansion are rigorously carried out within light-front QCD and without 
recourse to perturbation theory. The remaining problem is how to evaluate various matrix 
elements of bilocal currents. These matrix elements contain both hard and soft quark and 
gluon dynamics. As we have analyzed in this section, all the hard and soft dynamics probed 
through the structure functions are completely carried by the target's bound state in the 
present formulation. This is the main advantage of this formalism that allows us to explore 
the perturbative and nonperturbative contributions to the structure functions in the same 
framework. In the rest of this section, we shall propose a scheme for such an exploration. 

In Sec. V.A, the hadronic bound state is formally expressed in terms of Fock space 
expansion on the light-front by Eq. (|5.1|) , and it is determined in principle by the light-front 
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bound state equation of Eq. (|5.5| ). However, the difficulty in determining wave functions 



by solving Eq. (5^) is that the QCD Hamiltonian contains more than one energy scale. 
At different energy scales, QCD Hamiltonian can exhibit different aspects of the dynamics. 
Let us roughly divide the quark and gluon dynamics into two energy domains, namely, 
high energy and low energy. In the high energy domain, the dynamics is controlled by the 
renormalized QCD Hamiltonian with all the constituents carrying momenta greater than a 
scale /ifact (~ IGeV) which we call the factorization energy scale. This high energy QCD 
Hamiltonian describes all the hard dynamics of quarks and gluons and determines the hard 
contributions to the structure functions which can be calculated in the perturbation theory. 
In the low energy domain, the effective QCD Hamiltonian is still unknown but such a low 
energy QCD Hamiltonian should fairly determine the low energy structure of the hadrons 
and is responsible for the soft contributions to the structure functions. 



Schematically, we may write the QCD Hamiltonian on the light-front for DIS as 

Hq(jj~, = /fc2^>^2 dkf(Pki±HQQj^{ki) for hard contributions , 

Hq^j-, = J dkf d?ki^l-LQ(jj^{ki) for mixed hard and soft mod 

Hq(j^ = /fc2^<^2 dkfdPkixHQfjjjiki) for soft contributions , 



ttLF 
^QCD 



(5.19) 



where Hqq^ represents the canonical light-front QCD Hamiltonian (with density Ti-QQ^ 
given in 0) in which the transverse momenta of all the quarks and gluons are restricted to 
be < fc^ < (i.e., hard partons), and Hq^^j denotes a low energy effective light-front 
Hamiltonian in which all the constituents have the transverse momentum fc^ < nf^^^^ (soft 
partons). This low energy Hamiltonian is, in principle, obtained by integrating out all modes 
with fc^ > fif^^^ from the canonical light-front QCD Hamiltonian, which leads to 'Hqqj^. In 
addition, we also introduce a Hamiltonian Hq^j^ which depends only on the interaction 
part and which mixes the hard and soft partons. Writing the light-front QCD Hamiltonian 
in such three parts will make the discussion of the perturbative and nonperturbative QCD 
contributions to DIS structure functions much more transparent, as we will see next. 
Now, the target bound state can be expressed by 

|P5) = f/,|P5,/iLt), (5.20) 

with 

rO 



Uh = T+ exp I - ^ dx^iH^cD + H^cd) | , (5-21) 
HI,cd\PS, /iLt) = \PS, /^Lt) • (5.22) 



In Eq. (|5.21| ), and if*^ contain the interaction parts only and the mixed Hamiltonian 
^QCD is active only in the extreme right of the time-ordered expansion. In other words, 
the hard and the soft dynamics in the bound states are determined separately by -ffgcD 
and Hq^jj but these two contributions are connected by Hq(j^ through the time-ordered 
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expansion of Eq. ( p.21| ) on the state \PS,^'f^^^) in Eq. (|5.20|) , where the soft dynamics, 
represented by \PS, fij^^^^) , must be solved nonperturbatively from Eq. ( ^.22| ), and the key 
point to solve Eq. ( p.22| ) is to find the low energy effective Hamiltonian Hq^j^. A practical 
procedure to find -f^gc-D light-front may be the use of similarity renormalization group 

approach plus a weak-coupling treatment developed recently 10-0 • Indeed, a major effort 
on the study of light-front QCD is underway at present to solve this problem ||26[| . 

To see how the perturbative and nonperturbative QCD contributions can be separately 
evaluated in the present formalism and how these two contributions are connected by Hqqj-), 
we substitute Eqs. ( p.2(Jfj5.2^ ) into the expressions of structure functions. Denote the struc- 
ture functions simply by Fi =: {F^, F2, ^fi, Sfr}, 



(5.23) 



where Fj involves the Dirac 7-matrices [see Eqs. ( [4.12| )-( ^.18 )]. It follows that 



711,712 



X (ni|t/^-^[^jr)riVa(0) ± h.c]Uh\n2) , 



(5.24) 



where |ni), \n2) are a complete set of quark and gluon Fock states with momentum kf < 
nfg.cf This is indeed the generalized factorization theorem in the light-front Hamiltonian 
formulation. The hard contribution is described by the matrix element, 

{ni\U^'[^^{r)TiMO) ± h.c]Uh\n2) , (5.25) 

which can be evaluated in the light-front time-ordered perturbation theory |^ . The physical 
picture corresponds to the multi-parton forward scattering amplitude with all the internal 
partons carrying a momentum with the transverse component k±: fif^^^^ k\ < and the 
longitudinal momentum fraction y: x < y < 1. The soft contribution is characterized by 
the overlap of the multi-parton wave functions in different Fock states: 



(5.26) 



which contains all the quantum correlations and interference effects of multi-parton (quarks 
and gluons) dynamics in the low energy domain with k^ < yUfa^f Since all the internal 
partons in the time-ordered expansion of Uh in Eq. ( [5.25[ ) carry momenta fif^^^ < k\ < Q"^, 
the mixed Hamiltonian Hqqj-, has the contribution only in the extreme left and extreme 
right of the time-ordered products. It is this effect that connects the hard contribution of 
Eq. ( |5.25|) to the soft contribution Eq. ( ^.26| ). We will present more detailed discussion in 

i- 

The simple parton picture in deep inelastic processes corresponds to the case of |?t,i) = 
\n2) in Eq. ( ^.24| ) with only one parton in |ni) actively participating in the high energy 
process, all others being spectators. This immediately leads to 



1 

dy'Ppp' ^i{y , X, 2 

A'-fact 



)<?m(y,Ais 



act,' 



(5.27) 
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where the hard scattering coefficient Vpp'^i is determined by 



act 



Here we have denoted \y,k_i_,s) {y = / P^) as the active parton state. Eq. ( ^.28| ) means 



that we have suppressed all references to the spectators in the states |?t,i). The hard scat- 
tering coefficient is directly related to the so-called the splitting function whose physical 
interpretation is the probability to find a daughter parton p' in the active parent parton p. 
The quantity Qaiiy , fJ^fact) y usually called the parton distribution function, is given by 

qUy,^Lt)=T.\{PS,^Lt\n)\\ (5.29) 

n 

where n runs over all the Fock states containing the active parton with momentum fraction y. 
Theoretically, the parton distributions are determined by solving Eq. ( |5.22| ). Physically, they 
contain only the quantum correlations of multi-parton dynamics but no quantum interference 
effects. Example of such distribution functions is given by Eqs. ( |5.7| - |5.8| ) for F2{x)/x which 
manifestly exhibits the simple parton picture. For detailed calculations also see Ref. 

The above discussions indeed constitute a presentation of factorization scheme in the 
light-front Hamiltonian formulation. The leading hard contributions to the structure func- 
tions are given by the the hard scattering coefficient Ppp',i(2/, x, ) and a complete calcu- 

' ^fact 

lation of Vpp'^i based on the light-front time-ordered perturbative expansion of the multi- 
parton wave functions will be presented in a subsequent paper |Q. The evaluation of soft 
contribution to the structure functions, given by qai{x, fi'^^^^) remains for future investiga- 
tions of nonperturbative light-front QCD approaches to the hadronic bound states. Other 
higher order contributions can also be systematically evaluated from Eqs. ( 5.25| ) and (|5.26|) 



Thus, a unified treatment of both perturbative and nonperturbative aspects of deep 
inelastic structure functions in the same framework may emerge which permits one to over- 
come the obstacles in dealing with the nonperturbative QCD dynamics in OPE and field 
theoretical parton model approaches. 



VI. PHYSICAL INTERPRETATION OF THE STRUCTURE FUNCTIONS FROM 

SUM RULES 



In this last section, we shall explore the physical meaning of the deep inelastic struc- 
ture functions in our framework of light-front QCD. The physical meaning of the structure 
functions can be easily understood from the sum rules they obey. Some of them have been 
known for long time but others are new. Sum rules generally arise from the existence of 
conservation laws. First we consider the case of sum rules in unpolarized deep inelastic 
scattering for which a detailed consideration of the energy-momentum density in QCD is 
necessary. 



A. Energy-momentum tensor in QCD 



The symmetric, gauge-invariant energy-momentum tensor in QCD is given by 
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(6.1) 



The last term vanishes using the equation of motion. Formally, we spht the energy momen- 
tum tensor into a "fermionic" part 9^'^ and a "gauge bosonic" part 9^'^: 



2 



and 



9^''^-F^"'F\, + -g^^{Fx.,f, 



(6.2) 



(6.3) 



with F" = d^Axa — dxA'^ + gfabcAlAxc- To be consistent with the study of deep inelastic 
structure function which is formulated in — gauge, we shall work in the same gauge. 
We have, for the fermionic part of the longitudinal momentum density. 



9 



I++ 



(6.4) 



= -F+^F+^ = d+A'd+A' 



(6.5) 



Thus 



9++ = iipr^+d^ip + d+A^d+A\ (6.6) 
free of interactions at the operator level itself. The longitudinal momentum operator 



dx d x±9 



2^ . Q++ 



Next consider the transverse momentum density 

1 



with 



Ot-i = ^^^7+^V and 9+1^ = ^,ptfd+ip. 



For the Hamiltonian density, the fermionic part is given by 

0+- = ^+-(1) + ^+-(2) 

with 
and 



(6.7) 



(6.8) 



(6.9) 



(6.10) 



(6.11) 
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^+-(2) = itjj-^d+ij-. (6.12) 
Using the Dirac equation for the fermion, we find that 0'^^^^^ = 9^^'^'^\ Thus we have. 



9^-=ii:-f-d^ij = 2iij-^d^i:- (6.13) 
= 2^+^[a±.{id± + gA±) + j°m]^[a±.{id± + gA^) + -f^m\ip+ . (6.14) 

The gauge boson part of the Hamiltonian density is more comphcated [0]: 

= (d^Ai)' + 2gr"^AlAid'Ai + ^ AlAiAi^Ai 
+ 2gd^A\ {r'^^d^^c + 2(V^+)tT>+) 

+ / (^) {r'^^bd^K + 2(^+)tT>+) (ip) {r'^'A^d+A^, + 2(^+)tT>+) 

(6.15) 

where we have used the equation of constraint for the gauge field. 

Next, we discuss the physical interpretation of the deep inelastic structure functions on 
the basis of the sum rules they obey. 

B. Longitudinal momentum sum rule 

The content of the momentum sum rule is known for long time. For completeness, we 
shall rederive it in our framework. The sum rule is simply that if we add up the longitudinal 
momentum fractions carried by all the quarks, antiquarks, and the gluons (alternatively by 
the valence quarks, sea quarks and the gluons) in the nucleon we should get one. /^From the 
expression of F2 in terms of the plus component of the bilocal current matrix element given 
in Eq.( [4.12| ) we have 

[dxF,{x) = T.^l(P I I P) ■ (6-16) 



Formally, we can define the "gluon structure function" |jTO| 

F,^{x) = ^ J dve-^^^P I F+-(r)F+."(0) I P) , (6.17) 

so that, 

I' dxF^ix) = (^(^) (P I I P) ■ (6.18) 
Only if we assume = 1, one can obtain the momentum sum rule 
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[' dx[F2 + F^] = 1. (6.19) 
Jo 



Similarly, from Eq.( [4.13[ ) in terms of the transverse component of bilocal current matrix 
element, we have 

/ dxF,{x,Q') = ^^i_^e^(PS|^(0)7+(^ 5)V^.(0)|P^) 

^ j:el{PS\eii,\PS), (6.20) 



and the momentum sum rule can also be written as 

^ -{P I I P) = T775T75r(^ I \(^P + 1 I ^) = 1 • (6.21) 



2P+P] ' ' ' ' 2P+P] 



The sum rule given in Eg. (|6.16| ) means that F2 measures the longitudinal momentum 
distribution of quarks inside the hadrons, as known long time ago. From Eqs. (|6.16 ) and 



( |6.20| ) we observe that the hadron expectation value of the longitudinal and transverse 
momentum densities gives the same information, namely, the total longitudinal momentum 
fraction carried by the partons. 

We note here an apparent paradox that results when one ignores the essential complexi- 
ties carried by the state. The operators corresponding to the transverse momentum density 
explicitly depend on the interaction since = 5* — igA^ and 6^12 depends on ip' which in 
turn depends explicitly on the interaction. Since we know that is a kinematical operator 
this appears puzzling at first sight. Thus we expect the apparent dependence of P' on the 
interaction to be spurious. However this cannot be demonstrated at the level of operators 
alone. But this is not a serious problem since what really matters are the matrix elements. 

Indeed, our demonstration []T2| that the matrix element of the transverse component of 



the vector bilocal has the same parton interpretation as that of the plus component and hence 
the apparent interaction dependence in the former is completely spurious in turns directly 
tells us that the interaction dependence of the operator 9^12 is completely spurious. In that 
case an explicit evaluation of off-diagonal matrix elements in the Fock space expansion of 
states is involved. Similarly, an explicit demonstration shows that O^li has no interaction 
dependence at the level of matrix element, namely 



9+\ = Uiji-f+d'ip + g{ij+yA'ilj+, (6.22) 



but the matrix elements of the second term vanishes. Then, at the level of matrix elements 

e+' = iip-f'd+ilj. (6.23) 

This demonstration clearly shows that drawing conclusions by looking at the operator struc- 
ture is quite misleading in the case of operators that are twist three in the conventional 
definition. 

Since F2 involves quark charges in specific combinations, it does not give the direct test 
of the above momentum sum rule. To test the sum rule experimentally, one can combine 
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the data for both the electron-proton and electron-neutron deep-inelastic scatterings and 
assume that the sea is flavor symmetric, then 



dx 



1 



9(P+ 
5 1 



9P+P] 



J:{ps\oU\ps) 

a 

■j:{ps\ota\ps) . 



(6.24) 



This shows that ^ J dx (x) + F2"' (x) is the total longitudinal momentum fraction carried 
by the all the quarks in proton and neutron. If the quarks carry all the momentum, then 
we expect that 



dx 



F?{x) + Fr 



X] 



5 
9 



(6.25) 



Experimental data shows that the above integral is 0.28. In other words, as is well-known, 
there are half of the momentum in hadrons are carried by gluons or the see quarks if the sea 
is not flavor symmetric. 

C. Sum rule for gi 



Now we consider the sum rule for gi and its physical interpretation. Integrating gi over 
x, we simply have 



ri(g^ 



dxgi{x, 



1 



2S+ 



J2el{PS\ijMl^l5M0)\PS) 



Note that 



(P^|^J0)7+75^.(0)|P5) = Aq^'mS' 



(6.26) 



(6.27) 



where Aq^^ is the distribution function of chirality carried by all the quarks in the longitu- 
dinally polarized target and renormalized in the gauge invariant scheme. The above form 
can also be directly obtained from Eq. (|4.3|). Then, we have 



(6.28) 



If one uses the chiral invariant renormalization scheme, the flrst moment of gi also exhibits 
the anomaly contribution pl|] . We will discuss how this property manifests in the light-front 



Hamiltonian formulation in a separate publication. 

It is clear now that gi describes the distribution of chirality carried by the quarks inside 
the target (proton or neutron). Note the flrst moment Fi is usually called the proton's spin 
structure function. But one must also be aware that on the light-front, the plus component 
of axial current is the same as the third component of the quark helicity operator density on 
the light-front. Therefore, its expectation value is the same as the third component of the 
spin on the light-front (i.e., light-front helicity), which is not the same as the 2;-component 
of intrinsic spin deflned in the rest frame of the equal-time coordinates. On the light-front. 
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as is well-known, there is a very complicated relation between the light-front helicity and the 
intrinsic spin in the rest frame. This relation depends on the interactions in the fundamental 
theory. At present, one only knows the exact relation for free theory [Q. In other words, 
gi does not really measure the spin of proton. Simply calling Eq. ( |6.28|) a spin sum rule is 
misleading. 



D. Light-front helicity sum rule 

For the fermions, the intrinsic light-front helicity distribution function is given by 

1 



Aq{x,Q' 



87^5^ 



dT]e-'''^PS\ ^{r)r^ HO) + h.c \PS) 



(6.29) 



where = ^7^7^. This is the same as the chirality distribution function gi. 
We define the orbital helicity distribution for the fermion 



1 



dr]e~'"'^{PS I ^{rh^i{x^d^ - x^d^)^{0) + h.c. \ PS). (6.30) 



47CP+ 

For the gluon, the intrinsic light-front helicity distribution is defined as 



Ag{x,Q') 



ATT{P+yX 



q(0) 



PS). 



The dual tensor 



F^"" = -e'^^'P^F,, with = 2. 

2 ^ 

We also define the light-front orbital helicity distribution for the gluon as 

1 



(6.31) 



(6.32) 



A(7i(x,g^ 



47^P^ 



dr]e~'"^PS I [x^F+°(r)5M«(0) 

-x2F+"(r)5^^a(0)l I PS). 



(6.33) 



Note that all the above distribution functions are defined in the light-front gauge A'^ = 0. 
The light-front helicity operator is given by 



= \l dx-d^x' 



(6.34) 



where 9^^ is the symmetric energy momentum tensor. Explicitly, the fermion orbital helicity 
operator 



-^2 I „,,+t 



Jq{o) = i J dx d X ip 
and the fermion intrinsic helicity operator 



x'd^ - x'd' 



1 r 

j3.^ = -J dx'd^x^^^^Ti^^^. 



(6.35) 



(6.36) 
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The gluon orbital helicity operator 
1 



x\d+A^d^A^ + d+A^d^A^) - x\d+A^d^A^ + d+A^d^A^) (6.37) 



and the gluon intrinsic helicity operator 

Jlii) = \l dx-Sx' 



A^d+A^ - A^d+A^ 



The helicity sum rule for the nucleon target implies 
1 



1 

2' 



(6.38) 



(6.39) 



where A/" = 2(27r)^P+5'^(0). Thus we arrive at the sum rule obeyed by the helicity distribu- 
tion functions 



dx 



Aq{x, Q^) + Aqiix, Q^) + A^(x, + A^l(x, 



1 



(6.40) 



as a result of light-front helicity conservation. 



E. Sum rule for gx and the Burkhardt-Cottingham sum rule 

From Eq. ( |4.15| ) it follows that 

rl 1 /•+00 r 1 , 

/ dxg^{x, Q^) = — dx dr]e~^- ^ el(A,^ + -P+rA2a). 

JO lOTT J-oo J Q ^ 2 ^ 

Provided the bilocal form factor A2a does not have pathological behavior as 
have, 

1 i_ 

dxgi{x,Q'^) = -Aia{0). 

o 



Also we have, from Eq. ( |4.17| ), 



dxgT{x,Q^) = -Aia{0). 



Since, Qt = 9i + 92, it follows from Eqs. (|6.42|) and (|6.43|) that 

/ dxg2{x, Q^) = 
Jo 



(6.41) 

0, we 

(6.42) 
(6.43) 
(6.44) 



which is the Burkhardt-Cottingham sum rule. 

Recently, in the literature, there have been discussions about the validity of BC sum rule 
in perturbative QCD ||2^. Here we have shown the validity of the sum rule exactly up to 
the leading contribution in the l/q~ expansion without recourse to perturbation theory. 
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Obviously, the BC sum rule does not provide us any intuition about the physical picture 
of g2. Indeed, as we have pointed out the physical picture for g2 is not clear, since exper- 
imentally one directly measures gi and gx when the target is polarized in the longitudinal 
and transverse directions, respectively. The transverse polarized structure function is gr 
rather than g2. Eqs. ( |6.42| ) and ( |6.43| ) indicates that by averaging over x, the longitudinal 
and the transverse structure functions give the same result. This can be again regarded as 
a consequence of rotational symmetry. However, this does not implies that gi{x) and grix) 
are the same. 

To see clearly the intrinsic physical picture of (/^(a;, Q^), let us consider the target state 
being transversely polarized in the x-direction without loss of generality. Then we can simply 
express \PS) as a combination of the helicity up and down states: \PS^) = -^(l-P T)±|-P i)) 
for = ±M. It is easy to show that gx measures the helicity flip processes on the light-front 



A 

+ h.c\P-\). (6.45) 



Not that the quantity g2 is purely introduced in the Lorentz decomposition of the hadronic 
tensor W'^'^ for historical reasons and has no clear physical interpretation. Only gi and gr 
have the clear physical picture: gi measures the parton helicity distribution and gr measures 
the parton helicity flip effect, which is equivalent to the measurement of the effects of chiral 
symmetry breaking and therefore it involves more complicated intrinsic dynamics of quarks 
and gluons. A possible relation between gT and dynamical chiral symmetry breaking is 
explored in pO|], and more detailed theoretical and experimental investigations remain to be 



further carried out. 



F. Sum rule for twist four part of Fl 



^From Eqs. (|4l^ ) and (|4l^) it follows that F2{-x) = F2{x) and Fl=^{-x) = -Fl=^{x) 
Consider the integral 



X Jo X 

P+ f+OO 



[ ^ xdx [ dve-'^^Y.el{P l^ainh- 

J-oo ''a 

^)^^iO) - h.c. \ P)^ (6.46) 



where in the first equality, the symmetry property of Fl has been used. Interchanging the 



orders of x and y integrations and carrying out the integrations explicitly, we arrive at |25 



/^xS!^ = A^e^[(P I v;jO).(7"9+ - p^7-^9+)^.(0) I P)]. (6.47) 
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Identifying iip^'d^ip = Oq~ , the fermionic part of the hght-front QCD Hamiltonian density 
and itp'j'^ ip = 0'^^, the fermionic part of the hght-front QCD longitudinal momentum 
density, (see Eqs. ( |6.4| ) and ( |6.13| ) above), we arrive at the interesting relation: 



X 



2 



(0) 



\2 qa 



(0) I P) 



(6.48) 



We have used the fact that the physical structure function vanishes for a; > 1. We observe 

that the integral of — ^ is related to the hadron matrix element of the fermionic parts of the 
light-front Hamiltonian density. The above relation makes manifest the non-perturbative 
nature of the twist-four part of the longitudinal structure function. 

The fermionic operator matrix elements appearing in Eq. ( 6.48|) changes with as a 
result of the mixing of quark and gluon operators in QCD under renormalization. Next we 
analyze this problem of operator mixing and derive a new sum rule at the twist four level 
arising as a result of the conservation of energy-momentum tensor. 

We define the twist four longitudinal gluon structure function 



i(<y)' 



1 xP+ 



d^ e 



{P I (-)F+^-(r)F,- (0) + -g^-F^'^^nFx.aiO) \ P) 



Then if we assume Cq, = 1, we have 
dx 



{p I F+^^(r)PAt(o) I p) 



(6.49) 



X 



Q 



^ \p I r-(o) I P) - I ^^^(0) I P) 



(6.50) 



But, 



(P I ^+-(0) I P) = 2P+P- = 2{M^ + (P^)2) and (P | ^++(0) | P) = 2(P^ 



(6.51) 



where M is the invariant mass of the hadron. Thus we arrive at a new sum rule for the twist 
four part of the longitudinal structure function 



'r=4 



(6.52) 



VII. CONCLUSION 

In this paper, beginning with an inverse power expansion of the light-front energy of the 
probe in the framework of light-front QCD, we have arrived at the most general expression for 
the leading contributions to deep inelastic structure functions as the Fourier transform of the 
matrix element of different components of bifocal vector and axial vector currents. Although 
some of the expressions are already known, others are either completely new, such as the 
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expression for F2 in terms of the transverse component of bilocal current matrix element 
and the expression for F^, or generahzations of earher results in some specific Lorentz frame 
to an arbitrary Lorentz frame, e.g. the expressions for (or (72)- 

We have also provided a consistent and Lorentz invariant definition of twist based on 
light-front power counting. Using the light-front power counting, we found that quark-gluon 
coupling operators contributing to g2 are indeed twist-two when we look at their Lorentz 
invariant matrix elements. Only Fl has the leading twist-four contribution. No twist-three 
Lorentz invariant matrix element emerges in the leading contributions of the inclusive deep 
inelastic lepton-nucleon scatterings. The light-front power counting ^ naturally separates 
the space-time dimensions into two different scales measured in deep inelastic processes, the 
transverse scale and the longitudinal scale. Only the transverse scale has the mass dimension 
which determines the l/Q"^ suppressions. Therefore the light-front mass dimension seems to 
be more meaningful than the concept of twist in the analysis of the l/Q"^ suppressions. Light- 
front power counting analysis may eliminate some confusions in the discussion of higher- twist 
contributions to structure functions, especially, the transverse polarized structure function. 

We have also derived new sum rules for qt and Fl, which provide the physical picture of 
these structure functions. An important feature of the present formulation of deep inelastic 
processes is the fact that we have unified the treatment of soft and hard contributions to the 
structure functions in terms of multi-parton wave functions. The hard contributions can be 
easily calculated from light-front time-ordered perturbative expansion in high energy QCD 
IHJ^Q, while the soft contributions can be evaluated by the multi-parton wave functions 
below the factorization scale (~ 1 GeV), by solving the light-front bound state equation 
based on the recently developed nonperturbative renormalization group approach 0-0] or 
other approaches on light-front QCD . Further investigations along this direction are in 
progress and will be published in forthcoming papers. 
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